The quasiparticle spectral function in doped graphene 
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We calculate the real and imaginary electron self-energy as well as the quasiparticle spectral 
function in doped graphene taking into account electron-electron interaction in the leading order 
dynamically screened Coulomb coupling. Our theory provides the basis for calculating all one- 
electron properties of extrinsic graphene. Comparison with existing ARPES measurements shows 
broad qualitative agreement between theory and experiment. We also calculate the renormalized 
graphene momentum distribution function, finding a typical Fermi liquid discontinuity ai kp- We 
also provide a critical discussion of the relevant many body approximations (e.g. RPA) for graphene. 
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The great deal of current activity [ij in two dimensional 
(2D) graphene arises from its possible technological sig- 
nificance as a new 2D electronic material where carrier 
density can be controlled by an external gate voltage and 
from its fundamental significance as a novel 2D zero band 
gap semiconductor system with chiral linear Dirac-like 
electron-hole band energy dispersion. In particular, the 
chiral linear energy dispersion with the conduction and 
the valence band crossing at the "Dirac point" has led 
naturally to the interesting analogy with QED whereas 
the gate-induced tunability of graphene carrier density 
brings up the tantalizing analogy with Si MOSFETs. An 
important question in this context is the extent to which 
the Coulomb interaction between carriers will modify or 
renormalize the chiral linear electron-hole band disper- 
sion in graphene. 

In this Letter we theoretically consider carrier inter- 
action effects in extrinsic graphene by calculating the 
many-body self-energy, the quasiparticle spectral func- 
tion, and the renormalized momentum distribution func- 
tion of graphene in the presence of free carriers (i.e. for 
doped or gated graphene where carriers fill the 2D band 
upto a Fermi level Ep). Our results show that although 
the quantitative renormalization effects of interactions on 
the graphene single-particle properties are substantial, 
extrinsic graphene remains an effective 2D Fermi liquid 
"metal" , qualitatively preserving its non-interacting chi- 
ral linear band dispersion even in the presence of mutual 
Coulomb interaction. 

The quasiparticle spectral function, A(k,w), is a cen- 
tral quantity in the many-body physics of interacting sys- 
tems with A(k, Lj) = —2lmG(k,Lu) where G'(k, w) is the 
single particle (retarded) Green function for momentum 
k and energy uj (we use h = 1 throughout this paper). 
For the non-interacting bare system we immediately get 
AoihjUj) — 2ttS{u! — Esk + Ep) where Ssk — s"fk with 
k = |k| is the bare graphene linear band dispersion with 
s — ±1 denoting the conduction (-1-1) and the valence 
(—1) band, and 7 « 10^ cm/s is the band velocity. We 
will assume that the system is n-doped with electrons 



filling the graphene conduction band upto a free carrier 
density (n) dependent chemical potential or Fermi level 
given by Ep = jkp, where the Fermi momentum kp — 
(7rn)^/^. We have taken into account the spin and the 
valley degeneracy of graphene in obtaining the Fermi mo- 
mentum. The noninteracting spectral function Ao(k,ti;) 
being a delta function signifies that the band electron at 
momentum k has all its spectral weight precisely at the 
energy Ek = jk, i.e. the noninteracting particle exists en- 
tirely at the energy jk for a given momentum k. In the 
presence of interaction effects, the many-body self-energy 
function E(k, w) modifies the single-particle Green func- 
tion: G~'^{\i,uj) = GQ^CkjLu) — E(k,a;), and the corre- 
sponding interacting or renormalized spectral function is 
given by 



A(k,u;) = 



2ImE(k,tj) 



[uj-£,k- ReE(k, uj)Y + [ImS](k, uj)Y 



(1) 



where E(k, w) = ReE(k, w) -|-iImE(k, w) is complex, and 
£,sk = £sk — Ep. In general, ^(k, w) could be a compli- 
cated function of k and w, and there is no guarantee that 
it will have a delta function peak defining a quasiparticle. 
We note that J |^A(k, w) = 1 is a sum rule, guarantee- 
ing that the electron at momentum k exists in the whole 
energy space, but it may exist completely incoherently 
spread out over the whole w-space without any coherent 
structure (i.e. a delta function at fc = kp), indicating a 
complete failure of the Fermi liquid picture. Thus, the 
Fermi liquid theory applies only when the renormalized 
spectral function A{k — kp,Lu) at the Fermi momentum 
has a delta function peak, i.e. A{kp,u!) can be written 
as A{kp,uj) — 2'kZ6{w — £,1) 4- Ain{uj), where Z is the 
so-called "renormalization factor" , f ^ denotes the renor- 
malized quasiparticle energy (measured from the chemi- 
cal potential) , and Aia is the incoherent background spec- 
tral function. If A{kp,u!) does not have any delta func- 
tion peak at all, then the system is a non-Fermi liquid. 

In Fig. 1 we show our calculated interacting quasipar- 
ticle spectral function for extrinsic graphene at a fixed 
carrier density n = 10^^ cm^^. The calculations are car- 
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FIG. 1: Calculated graphene spectral function for different 
wave vectors (a) fe = 0, (b) k — kp, (c) k — l.bkp without 
any disorder effects, (d) Spectral function as a function of 
wave vector for different energies (MDC) and (e) as a function 
of energy for different wave vectors (EDC). In (d) and (e) we 
include an impurity scattering rate of Q.5Ef as explained in 
Fig. 3(b). We only show the spectral features above the Dirac 
point. Here all energies are measured from the Fermi level, 
and the Dirac point is at lu — —Ep- 



ried out at T = 0, but thermal effects are unimportant 
here since T/Tp <C 1 even at room temperatures. The 
following salient features of graphene quasiparticle spec- 
tral function are notable in Fig.l: (1) There is a well- 
defined delta- function quasiparticle peak a,t kp with a 
rather substantial spectral weight of Z « 0.85, indicating 
that 15% of the bare spectral weight goes into incoher- 
ent background. (2) The quasiparticle spectral function 
for k ^ kp shows, in general, broadened peak structures 
indicating damped quasiparticles. (3) The generic broad- 
ened double-peak structure at k ^ kp indicates that away 
from the Fermi surface, the renormalized graphene spec- 



tra would have two distinct energies — the second peak, 
which has been well-studied in the literature in both 2D 

and 3D Q interacting electron systems, is often re- 
ferred to as the "plasmaron" peak, indicating a coupled 
electron-plasmon composite excitation. (4) All spectral 
functions for k ^ hp have finite width corresponding to 
quasiparticle damping defined by the imaginary part of 
the many-body self-energy, ImE(fc, w). Note that in Figs. 

1 (d) and (e) disorder contributes to combining the double 
peak structure into a very broadened single peak. 

Since Fig. 1 giving the interacting graphene quasipar- 
ticle spectral function is the central result being pre- 
sented in this paper, we first discuss the importance and 
the implications of our calculated spectral function be- 
fore describing the details of our theory and other re- 
sults. First, extrinsic (i.e. gated or doped) graphene 
is a Fermi liquid with a well-defined undamped quasi- 
particle at kp. Thus, the chiral Dirac-like linear disper- 
sion of graphene band structure does not lead to any 
anomalous non-Fermi liquid behavior. Second, extrin- 
sic graphene has well-defined, but damped, quasipar- 
ticle peaks for all momenta. One direct experimental 
probe of the quasiparticle spectral function is the tunnel- 
ing spectroscopy, which has been studied extensively in 
GaAs-based 2D systems [1], but has not yet been stud- 
ied in graphene. The measurement that comes closest 
to studying j4(k, lu) in graphene is the angle-resolved- 
photo-emission spectroscopy (ARPES) 0, Q. Unfortu- 
nately, there are problems in directly comparing our the- 
oretical spectral function of Fig. 1 with the experimental 
ARPES results. One problem is that electron-phonon 
interaction induced many-body renormalization |i7,] and 
disorder also contribute to the graphene spectral func- 
tion. In Figs. 1(d) and (e), we show the theoretical re- 
sults corresponding to the experimental ARPES spectra 
— in particular. Figs. 1(d) and 1(e) respectively cor- 
respond to the so called MDC (momentum distribution 
curves) and EDC (energy distribution curves) of ARPES 
data. Since the actual theoretical spectral function. Fig. 
1, has very strong momentum (fc) and energy (w) de- 
pendence, a direct comparison with ARPES data will 
necessarily involve detailed instrumental issues involv- 
ing resolution, the {k,u;) regime of averaging, and in- 
strumental errors which are all well beyond the scope of 
the current theoretical work. A very cursory comparison 
between our results and the experimental ARPES data 
[1,0] show reasonable qualitative agreement, but short of 
large-scale data-fitting, one cannot make definitive quan- 
titative statements. One important point to note here is 
that the plasmaron structure (i.e. the additional peak) 
does not really show up in the MDC and EDC spec- 
tra since they carry small spectral weight compared with 
the main quasiparticle peak, particularly in the energy 
regime near Ep as observed experimentally 

We now describe the theory leading to our calcula- 
tion of the interacting quasiparticle spectral function de- 
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picted in Fig. 1. The self-energy, S(k,w), defining 
the spectral function through Eq. (1) is given, in the 
leading order dynamically screened Coulomb interaction 
approximation [3j 

Es(k, ia;„) = -/csT^ ^ Gq^s' (k + q, za;„ + ii^„) 



e{q, ivr, 



-F,,,(k,k + q), (2) 



where Vc{q) = jnq is the 2D Coulomb interaction 
with the background lattice dielectric constant k and 
Fss/(k, k') = (1 + ss' cos 6*^15;' )/2 is a sublattice over- 
lap matrix element arising from graphene band struc- 
ture, and e(g,a;) = 1-1- V"c(g)n(g,a;) is the dynamical 
RPA dielectric function for graphene with the irreducible 
electron- hole polarizability n((7, w), which has recently 
been calculated 0. We note that the irreducible self- 
energy approximation used in our theory corresponds to 
keeping the Coulomb interaction and the infinite series of 
electron-hole polarizability bubble diagrams in the self- 
energy calculation. 

An important technical point in the calculation of the 
self-energy is that, in general, two distinct types of field- 
theoretical divergence may appear in the theory: infrared 
(small momentum) and ultraviolet (large momentum). 
The infrared divergence arises from the 1/g long-range 
{q — *■ 0) divergence of the Coulomb interaction, and is 
regularized by our RPA theory through screening, i.e. 
the fact that T4(g)/e(g,ti') does not have a q in- 
frared divergence. While the infrared self-energy diver- 
gence, which is regularized by RPA screening, is generic 
to all problems involving long-range Coulomb interac- 
tion, the ultraviolet divergence, which arises from the pe- 
culiar graphene band dispersion, is specific to graphene 
in the sense that it does not occur in the corresponding 
paraboHc band 2D (or 3D) system [1,0]. The ultraviolet 
divergence in the graphene self-energy is fixed by real- 
izing that the linear Dirac dispersion of graphene only 
applies upto momenta of the order of inverse lattice con- 
stant, and therefore all momentum integrals should have 
an upper cut off fee ^ 1/^1 where a is the graphene lat- 
tice constant. We emphasize that there is nothing ad hoc 
or mysterious about this ultraviolet regularization since 
the band dispersion in graphene would deviate strongly 
from linearity for k > kc- We also note that the in- 
frared divergence associated with the q ^ singular- 
ity of the long-range Coulomb interaction and ultravio- 
let divergence associated with the q > kc ^ lattice 
wave vector scales have completely different origins, and 
while RPA screening fixes the infra-red divergence, the 
momentum cut-off fixes the ultraviolet divergence. Both 
infra-red and ultra-violet regularizations are important 
for graphene in contrast to 2D (or 3D) parabolic disper- 
sion interacting electron systems where only the infra-red 
divergence is germane. 
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FIG. 2: Calculated graphene self energy for k = l.Skp as a 
function of energy, (a) ReE and (b) |ImE| of correlation part, 
(c) Plasmon contribution to the imaginary part of self energy 
for several wave vectors. The Dirac point is at a; = —Ef- 



In Fig. 2 we show our calculated graphene self-energy 
for n = 10^^ cm~^ at fc = l.bkp- We note that the peak 
in ImS correlates with the dip in the ReS, and most of 
the spectral weight resides near dip of ReS (or equiva- 
lently, at the peak of ImS). We emphasize that, as is 
apparent from the detailed spectral function in Fig. 1, 
the various structures in the self-energy lead to specific 
features in A(k, u) — in particular, ReS and ImS control 
respectively the energy renormalization and the broaden- 
ing (or damping) of the quasiparticle. An important issue 
in this context is the precise contribution of the collective 
plasmon mode to ImS, which we have evaluated explic- 
itly as shown in Fig. 2(c). Only far away from the Fermi 
surface the plasmon contributions are strong, giving rise 
to the strong plasmaron peak in the spectral function. 

In Fig. 3 we show our calculated renormalized 
quasiparticle momentum distribution function n{k) = 

/-M i^^e^' ^) extrinsic graphene and disorder effects 
on the spectral function. The effect of electron-electron 
interaction on n{k) is rather obvious in Fig. 3(a). We 
note that the interacting momentum distribution func- 
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FIG. 3: (a) Renormalized momentum distribution function 
as a function of momentum for = 0.8 (solid line) and 0.2 
(dashed line), (b) Disorder effects on the fc = spectral func- 
tion for different impurity scattering rates of 0.1, 0.3, 0.5-Ef. 



tion has a discontinuity of relative size Z at k — k p , 
clearly establishing the Fermi liquid behavior of extrin- 
sic (i.e. doped or gated) graphene. The discontinuity 
increases as Ts decreases. We note here that the corre- 
sponding intrinsic (or undoped) graphene does not have 
a discontinuity in the distribution function since it is a 
marginal Fermi liquid 

Impurity effects are usually introduced diagrammat- 
ically into the RPA screening by including ladder im- 
purity diagrams into the electron-hole bubble. We use 
a particle-conserving expression [loj. which captures the 
essential physics of impurity scattering effects on the elec- 
tron polarizability 11. We find that disorder strongly sup- 
presses plasmaron peak, but not the quasiparticle peak. 
Thus the double peak structure in the spectral function 
becomes a very broadened single peak. In the strongly 
disordered experimental graphene samples therefore we 
only have a broadened quasiparticle peak. 

We now comment on the experimental implications, 
the theoretical approximations, and the connections to 
earlier work. As mentioned before, the graphene ARPES 
measurements have information on the electron spectral 
function, which are, however, complicated by the pres- 
ence of additional interaction effects such as electron- 
impurity (i.e. disorder) and electron-phonon interac- 
tions. If these two interaction effects can be subtracted 
out (e.g. ~ the 200 meV structure ^ presumably arising 
from phonons in the experimental data), then ARPES 
measurements could indeed be compared with our calcu- 
lated spectral function shown in Fig.l. We note that, in 
agreement with the ARPES data, the quasiparticle spec- 
tral peak becomes wider due to plasmon contribution as 
the energy approaches to the Dirac point (lo = —Ep). In- 
cluding disorder effects, indeed there is good agreement 
between theory and experiment. 

Our theory is based on the leading-order expan- 
sion of the electron self-energy using the dynamically 



screened Coulomb interaction (the so-called infinite bub- 
ble diagram expansion with each bubble being the non- 
interacting electron-hole polarizability [1]), which we be- 
lieve to be a quantitatively accurate approximation for 
extrinsic graphene (i.e. at any finite carrier density) 
by virtue of graphene having a reasonably small (and 
density independent) dimensionless interaction parame- 
ter Ts = e^/K7 ~ 0.8 (0.4) for graphene on Si02 (SiC) 
substrate. In fact, the same self-energy approximation, 
often referred to in the literature as "RPA" or "GW" 
approximation, is known to work well in 3D alkali met- 
als where ^ 3 — 6 Q and in 2D semiconductor sys- 
tem [3, 0] where rg ~ 1 — 10. For extrinsic graphene 
only the terms having bubble diagrams (or RPA type di- 
agrams) in self-energy have the infra-red divergence. All 
other vertex corrections (i.e. non-RPA type diagrams) 
have finite values. Thus, most contributions (for < 1) 
come from RPA type diagrams. Since we can control 
the ultra-violet divergence by introducing a physical en- 
ergy cut-off, our RPA approach is quantitatively accurate 
for extrinsic graphene. However, for intrinsic graphene, 
there is no infra-red divergence in RPA type diagrams, 
only the ultra-violet divergence. Since all diagrams then 
contribute to the self-energy, the perturbative expansion 
does not converge for intrinsic graphene. Therefore, RPA 
is not a meaningful approximation for intrinsic graphene, 
and indeed it predicts a non- Fermi liquid behavior 



In discussing connection between our work and earlier 
work, we mention that, although there has been a great 
deal of recent theoretical work on interaction effects in 
graphene, much of it has focused on intrinsic graphehe 
111, and most of the work on extrinsic graphene has fo- 
cused on thermodynamic properties [l^ . A recent calcu- 
lation of the quasiparticle spectral function in graphene 
due to electron-phonon interaction has appeared FH| in 
the literature, and within the weak-coupling theory (i.e. 
both electron-electron and electron-phonon interactions 
being weak), the total spectral function should be a sum 
of these two spectral functions. We do note that our 
calculated spectral function and self-energy in extrinsic 
graphene is qualitatively rather similar to that in ordi- 
nary parabolic band 2D carrier systems , thus pointing 
to the fact that, in spite of its chiral Dirac- like band dis- 
persion, doped graphene in the presence of free carriers is 
qualitatively similar to doped 2D semiconductor systems. 



In summary, we have calculated the electron spectral 
function in 2D doped (i.e. extrinsic) graphene, finding 
our theory to be in good qualitative agreement with the 
available experimental data. Our theory convincingly 
demonstrates the Fermi liquid nature of doped graphene. 
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